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Abstract 

Surface plasmon polaritons (SPPs) in a semi-bounded degenerate plasma (e.g., a metal) are 
studied using the quasiclassical mean-field kinetic model, taking into account the spatial dispersion 
of the plasma (due to quantum degeneracy of electrons) and electron-ion (electron-lattice, for 
metals) collisions. SPP dispersion and damping are obtained in both retarded (u)/k z ~ c) and 
non-retarded (u)/k z <C c) regions, as well as in between. It is shown that the plasma spatial 
dispersion significantly affects the properties of SPPs, especially at short wavelengths (less than 
the collisionless skin depth, A < c/uj pe ). Namely, the collisionless (Landau) damping of SPPs (due 
to spatial dispersion) is comparable to the purely collisional (Ohmic) damping (due to electron- 
lattice collisions) in a wide range of SPP wavelengths, e.g., from A ~ 20 nm to A ~ 0.8 nm for 
SPP in gold at T = 293 K, and from A ~ 400 nm to A ~ 0.7 nm for SPPs in gold at T = 100 K. 
The spatial dispersion is also shown to affect, in a qualitative way, the dispersion of SPPs at short 
wavelengths A < c/u pe . 

PACS numbers: 52.35.-g,73.20.Mf 



* Electronic address: [y. tyshetskiy ©physi cs . usyd . edu.au 



1 



I. INTRODUCTION 



As early as 1950s it has been demonstrated theoretically [lj, and later confirmed experi- 
mentally 2H5J that bounded metallic structures (e.g., thin films), as well as other bounded 
plasmas 6- J, can support a special type of collective oscillations of the plasma electrons, 
surface plasma waves, also called surface plasmons (SP), that propagate along the plasma 
boundaries and whose field and energy density are concentrated near the boundaries. 

Since then, there has been a significant advance in theoretical and experimental investiga- 
tions of surface plasma waves and their applications in various bounded plasma structures, 
both in the field of plasma science (see Refs js-ll] and references therein) and in the fields 
of condensed matter and surface science (see, e.g., Ref. 12j for a review). Currently, there 
is a renewed interest in surface plasmons due to their ability to concentrate light in sub- 
wavelength structures, enabling to create surface plasmon-based circuits that can couple 



photonics and electronics at nanoscale 



13H15|. This offers a route to faster and smaller de- 



vices, and even to new technologies employing surface plasmons 16j. For example, one of the 
recent interesting advents in the new area of quantum nanoplasmonics is the development 
of the concept of spaser (a surface plasmon "laser") [17|, followed by its further develop- 



ment into a lasing spaser 



nanolaser 
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and by an experimental demonstration of a spaser-based 



20). 



In view of these developments, understanding the properties of surface plasmons in various 
metallic (and semiconductor) structures, bounded by vacuum or dielectric, is thus important. 
This requires using models for the dynamic response of charge carriers in such structures 
to self-consistent electromagnetic fields, that appropriately take into account the relevant 
effects arising from quantum nature of the charge carriers in such structures and from their 
interaction with the underlying ion lattice. Such quantum effects may significantly alter the 
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221. 



properties of the surface waves; see, e.g., Refs. 

In general, the electric field of a surface plasmon has both longitudinal and transverse 
components, i.e., a surface plasmon is, in general, an electromagnetic wave coupled with a 
collective oscillation of surface charge; such a hybrid electromagnetic surface wave is called 
a surface plasmon polariton (SPP). The properties (dispersion and damping) of SPP in a 
semi-bounded degenerate plasma (e.g., a metal) are a subject of this paper. 

In the non-retarded limit, when the SPP phase velocity is much smaller than the speed 
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of light, the SPP field becomes purely electrostatic. In this limit, the SPP corresponds 
to a longitudinal surface charge density wave; its properties are thus affected only by the 
longitudinal part e of the plasma dielectric response. Recently, the properties of such purely 
electrostatic surface plasmons (ESP) in a semi-bounded plasma with degenerate electrons 
were analyzed using the semi-classical kinetic model 2^|. In particular, it was shown that 
the spatial dispersion of degenerate plasma has a significant effect on the frequency and 
damping rate of ESP at short wavelengths, k z Xp > 1, where k z = 2n/X is the ESP wave 
vector component along the plasma boundary, A is the ESP wavelength, and A^ = vf/V3i 



pe 



is the Thomas- Fermi length (vf = h\/3ii 2 n e /m e is the Fermi velocity of plasma electrons, 
co pe = (A-Ke 2 n e /m e ) 1 / 2 is the electron plasma frequency, — e and m e are the electron charge 
and mass, h is the reduced Planck constant, and n e is the electron number density in plasma; 
we use CGS units). 

In this paper, we aim to generalize the work on electrostatic surface plasmons in a semi- 



bounded collisionless degenerate plasma [23] , and consider the properties of electromagnetic 
surface plasmons, or surface plasmon polaritons, in a semi-bounded degenerate plasma with 
electron- ion collisions (e.g., SPPs in a semi-bounded metal where the electron-lattice colli- 
sions play a significant role). In particular, we consider effects of plasma spatial dispersion, 
and collisions, on SPP dispersion and damping. We show that the collisionless damping 
of SPPs (due to spatial dispersion of a metal) is comparable, or even exceeds the purely 
collisional (Ohmic) damping of SPPs (due to electron-lattice collisions) in a wide range of 
wavelengths, e.g., from A ~ 20 nm to A ~ 0.8 nm for SPPs in gold at T = 293 K, and 
from A ~ 400 nm to A ~ 0.7 nm for SPPs in gold at T = 100 K. The spatial dispersion 
is also shown to affect the SPP spectrum, especially at short wavelengths (less than the 
collisionless skin depth, A < c/u pe ), where the spatial dispersion changes SPP spectrum in 
a qualitative way. We conclude that the spatial dispersion (due to quantum degeneracy of 
plasma electrons) is rather important for damping, and, to a lesser extent, for dispersion of 
SPPs in a semi-bounded degenerate plasma (e.g., a bounded metal), and should be taken 
into account in the relevant SPP models. 
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II. METHOD 



A. Model and Assumptions 

We consider a semi-bounded, nonrelativistic plasma (e.g., a metal) with degenerate 
mobile electrons (T e <C e F , where T e is the electron temperature in energy units, e F = 
h 2 (3ir 2 n e ) 2 / 3 / 2m e is the electron Fermi energy), and immobile ions; the equilibrium number 
densities of electrons and ions are equal, no e = n-oi = no (quasineutrality). The plasma is 
assumed to be confined to a region x < 0, with mirror reflection of plasma particles at the 
boundary x = separating the plasma from a vacuum at x > 0. In the absence of fields, 
the equilibrium distribution function of plasma electrons /o e (p) is defined by an isotropic 
Fermi-Dirac distribution, which in the limit T e <C e F reduces to 

r p 2 /2m e - e F (n e )' 



(2tt^) 



-a\p F (n e ) -p] , (1) 



T 

where p F (n e ) = \/2m e eF(n e ) is the electron Fermi momentum, <j(x) is the Heaviside step 
function. 

Following the discussion of Ref . 23( , we adopt here the quasiclassical kinetic description 



of plasma electrons in terms of the 1-particle distribution function / e (r, p, t) 24j], whose 
evolution is described by the kinetic equation 

9 f< + JL . °fl _ e ( E + 1211 V ?1± = V IeMe: U (2) 



dt m P dv 



where I ea {fe, fa) describes collisions of electrons with particles of sort a = e,i, i.e., with 
electrons and ions. Here we are only interested in the electron-ion collisions (due to coherent 
scattering of electrons on the ions, which in metals is equivalent to electron-phonon scatter- 
ing), which transfer the energy of quiver motion of electrons in the field of collective plasma 
oscillations (in particular, in the field of SPP) to the "ionic thermal bath" (e.g., ion lattice 
in metals), resulting in SPP damping. 
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B. Electron-ion collisions 



In general, collisions between degenerate electrons and non-degenerate classical ions can 



be described by the quantum Lennard-Balescu collision integral 
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26] 



PJ 



y ei (p - p')r 5 k(p') + Q(pD - e e (p) - 

XS [ P + Pi - P ' - P J] {/e(p') [1 - /e(p)] /i(pj) - /e(p) [1 " /e(p')] /.(ft)} , 



(3) 



where e e ,i{p) — |p| 2 /2m ej j, and is the Fourier transform of the dynamically screened 

electron-ion interaction potential, 

4irh 2 e 2 Zi 

VM) = / ^exp (--q-rj V ei (r) = 



where e l (u,k) is the longitudinal dielectric response of the plasma, 
hr 1 [e e (p + q) ~ e e(p)] ; Q is the momentum transfered in the collision, and Zj is the ion- 
ization number of plasma ions. 

In a weak spatially uniform oscillating field, E = Eo exp(— iut), the electron distribu- 
tion function can be approximated as / e (r,p, t) = fo e (p) + 5f e (r,p,t), with a perturbation 
Sf e (r, p, t) = g(r,p) cos9, where 9 is the angle between p and E, and where \g(r,p)\ <C fo e (p) 
is assumed (due to the field being weak) . In this case the electron-ion collision integral ([3]) 



reduces to 



26] 



<*/ c (r,p,t) 



fe — f{ 



T ( \ - ^J_^jJLl ' - Je ~ JOe (A\ 
leiyP) — i \ I \ i W 

Teiip) T ei (p) 

where r~* (p) is the electron- ion collision frequency for electrons with momentum p, defined 



as 
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22] 



^(P) 



n 0i m e 



2p ' h dk 



a(k)k 3 



AnZ.e 2 



k 2 e l {0,k) 



(5) 



2p 3 J 2tt 

Note that ([5]) is obtained by assuming static screening, e(u,k) ~ e(0, k), for simplicity. 
The structure function a(k) in (J5j) depends on the arrangement of ions in the system con- 
sidered, and can be obtained from experimental measurements of x-ray scattering on the 



corresponding medium 



Yet in case of a weak non-uniform oscillating field (e.g., the self-consistent field of plasma 
collective oscillations, such as SPPs), the perturbation 5f e (r, p, t) is no longer proportional 
to cos#, and thus the integral fl3]) no longer reduces to the simple form (J4]). Nevertheless, 
the form (jl]) is rather attractive due to its simplicity, and it is thus desirable to use this 
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form in our model, especially since it becomes exact for fields that are close to uniform, 
i.e., for plasma perturbations with k\F <C 1. We therefore adopt a model collision integral 
describing collisions of electrons with the heavy plasma ions (e.g., with the ion lattice in 
metals) in the form of Bhatnagar-Gross-Krook (BGK) integral 29J: 



tBGK , 



P) = -f« [fe - foe(n e ,T e )} 



(6) 



where fo e (n e , T e ) is the quasi-equilibrium distribution of electrons given by <Q with ep(n e ) = 
h 2 (37r 2 n e ) 2 / 3 /2m e , in which n e and T e are the perturbed electron density and temperature, 
defined as 



n e (r,t) = J f e (r,p,t)d 3 p, 



T e (r,t) 



v 2 f e (r,p,t)d 3 p. 



(7) 
(8) 



2n e (r,t) 

The model collision integral (EJ) conserves the number of particles, and energy, and mimics 
the relaxation of electron distribution towards the quasi-equilibrium distribution fo e (n e ,T e ), 
occurring over the characteristic time . The value of v e { can be estimated as v e i ~ 
1/T ei (pp), with T e i(pp) from ([5]), for which the knowledge of the structure function a(k) is 
required. However, for a crude estimate (sufficient for our model BGK integral), we can 
simply use the dc resistivity of the bulk plasma, p = v e im e /e 2 rio e , which is available for most 
metals from experimental measurements. This gives 



Eouipep (in SI units). 



(9) 



pe 



For example, for gold at room temperature (T = 293 K) p ?a 2.2 x 10 8 Ohm-m 30|, 



u pe w 1.3 x 10 16 s" 



30 



3l{ , which gives 



pe 



2.5 x 10' 



(10) 



Au, T=293 K 



C. Dispersion equation for SPP 



Introducing a small perturbation 5f e — f e — fo e associated with a weak electromagnetic 
field E, B, and linearizing the kinetic equation (j2J) with I ee = and I e i defined by (Q, we 
obtain the equation for Sf e in the form 



A + . A ) 

dt m e dr 



eE- 



dfi 



0e 



dp 



Sfe + 



2e F (n ) df { 



0e 



3 n Q de 



no 



Sfed'p' 
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where we have used the isothermal approximation for the BGK collision integral (see Ap- 
pendix |A|). Eq. ffTTj) is complemented with the Maxwell's equations, with charge and current 
densities defined in terms of 5f e as p q = —e f 5f e d 3 p, j = — (e/m e ) f p5f e d 3 p, with im- 
posed continuity of E, B components parallel to the plasma boundary x = 0, across the 
plasma-vacuum interface. Assuming a specular reflection of plasma particles (electrons in 
our case) at x = 0, and following the standard procedure (see, e.g., Ref. 29]), we obtain the 
dispersion relation for TM-polarized SPP (note that TE-polarized SPPs are not allowed in 
the considered system): 



k 2 



+0 ° dk x 

.. — , rh 



ki 



2U2 



U 2 k 



0. 



(12) 



.e'(w,k) c 2 k 2 -u 2 e tr (co,k) 

where u and k are the frequency and wave vector of TM SPP, respectively, k x and k z are 
the components of k perpendicular and parallel to the boundary x — 0, respectively, so that 
k = |k| = ^k 2 x + k 2 . Note that in the non-retarded limit, u/k c, the dispersion equation 
(TI2"]) for SPP reduces to the dispersion equation for ESW, studied in Ref. [231 ] : 



k y 

1 + — 

7T 



0. 



, k 2 e l (u,k) 

The plasma properties with respect to the weak electromagnetic field of SPP are entirely 
contained in the longitudinal and transverse linear dielectric response functions e l (u, k) and 
e tr (u, k), which, for the model (1TTT) with isotropic equilibrium distribution ([T]), are found to 



be 



29] 



e l {u,k) 



3 p 



k 2 Vp L 



W + W e i i /" W + Wei + fcu F 

1 m 



2kvp 
u + iu ei + kvp 



n -1 



e ir (w,A;) 




(13) 



(14) 



III. PROPERTIES OF SURFACE PLASMON POLARITONS 



Eq. f|T2l with e l (u,is) and e tr (o;,k) defined by Eqs ffl~3"]) - fn4"|) describes dispersion and 
damping of SPP in a semi-bounded degenerate plasma (e.g., a metal) with electron-ion 
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collisions approximated by the BGK integral in the right-hand side of (llip . Below we analyze 
the effects of plasma spatial dispersion and electron-ion collisions on SPP properties. 

To see how the spatial dispersion (fc- dependence of plasma responses £ l (co, k) and e tr (u, k), 
due to electron velocity spread vp arising from Pauli blocking) and electron-ion collisions 
affect the dispersion and damping of SPPs in the considered system, it is instructive to 
compare the following three models: 

1. "Local collisional" approximation for the medium, in which the electron- ion collisions 
are retained, while the spatial dispersion is neglected, by formally taking the limit 
vp — > in fri3l) - ffT^l) . which gives 

,2 



£\u,k) = e tr (u,k) 



id. 



pe 



(15) 



lo (u + iu ei ) ' 

This model describes purely collisional effects on SPP dispersion and damping. The 
SPP frequency and damping rate obtained in the "local collisional" approximation are 
denoted by u VF=0tUei>0 and 7<, F =o,!^>o, respectively. 

2. "Nonlocal collisionless" approximation, in which the spatial dispersion is retained, 
while the electron-ion collisions are neglected by taking the limit v ei — > in ffT3]) - 
( fl4l) . This model describes the effects of spatial dispersion on SPP properties, e.g., 
collisionless damping of SPPs. The SPP frequency and damping rate obtained in 
the "nonlocal collisionless" approximation are denoted by u} VF >o,v ei =o and ■y VF >o,u ei =o, 
respectively. 

3. "Nonlocal collisional", i.e., the full model, in which both spatial dispersion and colli- 
sions are retained in (flBl -t CHl) . This model describes the combined effects of spatial 
dispersion and collisions on SPP properties. The SPP frequency and damping rate 
obtained in the "nonlocal collisional" approximation are denoted by w 1 , F >o,i/ ei >o an d 
lv F >o,u et >o, respectively. 



Only the weakly damped SPPs are of physical interest. T 
are obtained in the same way as described in detail in Ref. 23] 



reir dispersion and damping 
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FIG. 1: (Color online) SPP spectra for gold at room temperature (T = 293 K), obtained from the 
"local collisional" model (solid blue curve), "nonlocal collisional" (full) model (solid red curve), 
and "nonlocal collisionless" model (not seen due to exact matching with the red curve). The light 
dispersion ui = ck z is shown with the dotted black line, the "cold plasma limit" oj = u pe /\/2 is 
shown with the dotted blue line, and the zero sound asymptote (I17h is shown with the dotted 



30]: uj pe = 1.3 x 10 16 s" 1 , v ei /uj pe = 2.5 x 10~ 3 [see Eq. ([TUP ]. 



green line. The gold parameters are 
vf/c = 5 x 10~ 3 . 

A. Dispersion 

The spectra of SPP u(k z ) obtained from the three models introduced above ("local colli- 
sional" , "nonlocal collisionless" , and "nonlocal collisional" ) are shown in Fig. HJ for the case 
of gold at room temperature. It is seen that the SPP phase velocities given by all three 
models are smaller than the speed of light, u/k z < c; i.e., the surface plasmons in a semi- 
bounded plasma with an ideal boundary are essentially non-radiative (they cannot decay by 
emitting a photon, and cannot be excited by a photon incident on the plasma surface). The 
dispersion curves for the "nonlocal collisionless" and "nonlocal collisional" models match, 
hence the collisions do not affect the SPP spectrum. The effect of the spatial dispersion 
on SPP spectrum, however, is significant, especially at short SPP wavelengths (A < 15 nm 
for gold at room temperature), as seen from Fig. [TJ The asymptotic expressions for SPP 
frequencies u(k z ), given by the "local collisional" and "nonlocal collisional" models match: 



w df=o,^„>o — ^F>o,^i>o — ck z — O [(c/c z ) 2 ] , ck z 
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i.e., the spatial dispersion does not affect the SPP spectrum significantly at large wave- 
lengths. However, at short wavelengths, the "local collisional" and "nonlocal collisional" 
models yield qualitatively different asymptotes for SPP frequencies: without the spatial 
dispersion, the SPP frequency approaches the "cold plasma limit" 



Uv F =0,u ei >0 = Ck z > W pe , (16) 

while accounting for the spatial dispersion leads to a qualitatively different SPP frequency 



asymptote 



23|: 



^v F >0,v ei >0 



k z v p ( 1 + 2 exp 



4 k 2 z v 2 F 

3 Wpe 



k z vp, ck z 3> uj, 



pe • 



(17) 



It is interesting to note that this part of the SPP spectrum (fTTj) is exponentially close to the 
spectrum of volume zero sound mode, propagating along the boundary (with k = k z ) in an 
uncharged Fermi gas. 



B. Damping 

The SPP damping rates 7(& z ), obtained from the three models ("local collisional", "non- 
local collisionless" , and "nonlocal collisional") are shown in Fig. [2j for the case of gold at 
room temperature. The purely collisional (ohmic) SPP damping rate, given by the "local 
collisional" model, has the following asymptotes: 

7^=o^ ei >o = — = a - , ck z < u pe , (18) 

'y VF =o,u ei >o = ck z ^>u pe . (19) 

The spatial dispersion of the plasma (due to quantum degeneracy of electrons) leads to 
collisionless SPP damping, which exhibits a non-monotonic dependence on k z , as seen in 
Fig. |2j at small k z , the collisionless damping rate monotonically increases with k z , reaching 
a distinct maximum of •y VF >o,u E1 =o ~ 6.2 x 10~ 3 u) pe at ck z /u pe ~ 0.7c/vf, and then decreases 
monotonically with k z for ck z /u pe > 0.7c/ Vp, quickly approaching zero at ck z ^> u pe . 

It might seem somewhat counter-intuitive that SPPs are subject to collisionless damping 
at all, despite the fact that their phase velocity is larger than the maximum electron velocity 
in degenerate plasma, u/k z > Vp- Indeed, it is well known that the volume plasmons in 
degenerate plasma, described within the quasiclassical model employed here (with quantum 
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J ' ' ' ' Ck- /W pe 

0.01 0.1 1 10 100 

FIG. 2: (Color online) SPP damping for gold at room temperature (T = 293 K), obtained from 
the "local collisional" model (solid blue curve), "nonlocal collisionless" model (solid green curve), 
and "nonlocal collisional" (full) model (solid red curve). The dotted blue line marks the maximum 
purely collisional damping 7 = — f e j/2 [see Eq. (fT9|) ]. The parameters are the same as in Fig. [TJ 

recoil ignored), are not subject to collisionless damping at all, exactly due to their phase 
velocity being larger than the maximum electron velocity in degenerate plasma, u/k > vf, 
so that no plasma electrons can be in resonance with the plasmons. However, a surface 
plasmon is in fact a result of interference of many (a continuum, in case of SPP in a semi- 
bounded plasma) "virtual" volume plasmons with a given k z , but with all possible k x ranging 
from minus to plus infinity [see Eq. ( TT2]) ]. As a result, the surface plasmon consists of virtual 
volume plasmons with phase velocities that are both larger and smaller than vp , and the 
latter are subject to resonant collisionless damping on plasma electrons, resulting in a finite 
collisionless damping of the resultant surface plasmon, even despite its phase velocity u/k z 
exceeding the electron Fermi velocity. 

As seen from Fig. [2], for gold at room temperature (T = 293 K), as well as for other metals 
with similar properties (e.g., silver, aluminium), the collisionless damping rate is equal to, 
or larger than the purely collisional damping rate, for SPP wavelengths ranging (for the case 
of gold) from extreme ultraviolet (ck z /u pe ~ 6, corresponding to A ~ 20 nm) to soft x-ray 
(ck z /uj pe ~ 2 x 10 2 , corresponding to A ~ 0.8 nm). The full SPP damping rate / jv F >o,iy ei >Oi 
due to the combined effect of collisions and spatial dispersion, given by the full "nonlocal 
collisional" model, is approximately (for v ei /u pe <C 1, which is the case for metals) a sum of 
the purely collisional 7„ F =o,i/ ei >o an d the purely collisionless 7„ F >o, I / (i4 >o damping rates, and 
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ck. /a>„ 



FIG. 3: (Color online) SPP damping for gold at T = 100 K (p = 6.5 x 10~ 9 Ohm m [3fl|, 
corresponding to v e i/uj pe ~ 7.5 x 10 -4 ), obtained from the "local collisional" model (solid blue 
curve), "nonlocal collisionless" model (solid green curve), and "nonlocal collisional" (full) model 
(solid red curve). The dotted blue line marks the maximum purely collisional damping 7 = —v e i/2 
[see Eq. (Tl9|) ]. 

is dominated by the collisionless damping at short SPP wavelengths, as seen from Fig. |2j 
At lower metal temperatures, the resistivity and, correspondingly, v e i drop 30], resulting in 
widening of the wavelength range at which the SPP collisionless damping dominates over 
the collisional damping. For example, at T = 100 K, we have p = 0.65 x 10 -8 Ohm m 30], 
corresponding to u ei /u pe ~ 7.5 x 10~ 4 (see Eq. Q), and the range of SPP wavelengths 
at which the collisionless damping dominates extends all the way from A ~ 400 nm to 
A ~ 0.7 nm, as seen in Fig. [31 



IV. CONCLUSION 



In this paper, surface plasmon polaritons in a semi-bounded degenerate plasma (e.g., a 
metal) were studied, using the quasiclassical kinetic model with the BGK collision integral 
modeling electron-ion (electron-lattice for metals) collisions. Spectrum and damping rate 
are obtained for SPP in both retarded and non-retarded cases, thus extending the results 
of Ref. 23| in which only the non-retarded limit was treated. The effects of plasma spatial 



dispersion (due to velocity spread of degenerate electrons) and of electron-lattice collisions 
on SPP dispersion and damping are studied. It is found that the collisions do not affect SPP 
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spectrum, but contribute to SPP damping. The spatial dispersion, however, significantly 
affects both spectrum and damping of SPPs, especially at short wavelengths. Namely, it is 
due to the plasma spatial dispersion that the SPP spectrum changes from the "cold plasma 
limit" to = co pe / \/2 to the "zero sound" asymptote to ~ k z vp at short wavelengths, A < c/co pe . 
In metals, the damping of SPP is approximately a sum of purely collisional (Ohmic) and 
collisionless (Landau) damping, and is dominated by the spatial dispersion at a wide range 
of wavelengths, especially at low metal temperatures. For example, in gold at T = 293 K, 
the total SPP damping is dominated by the collisionless damping (due to spatial dispersion 
arising from quantum degeneracy of conduction electrons) for SPP wavelengths ranging from 
A ~ 20 nm to A ~ 0.8 nm, while in gold at T = 100 K (e.g., cooled in liquid nitrogen) the 
total SPP damping is dominated by the collisionless damping in a much wider range of SPP 
wavelengths, spanning from A ~ 400 nm to A ~ 0.7 nm. We thus conclude that the spatial 
dispersion is important for SPP properties, especially in metals at low temperatures, and 
should be taken into account in the relevant models. 
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Appendix A: BGK integral in isothermal approximation 

The BGK collision integral is defined by (EJ) with fo e (n e ,T e ) being the quasi-equilibrium 
Fermi-Dirac distribution depending on the perturbed electron density n e and temperature 
T e , which for T e <C epijie) is given by (Tj[|). At equilibrium we have f e = foe(n 0e ,T 0e ), 
where riQ e and T 0e are the equilibrium electron density and temperature. Introducing a 
small perturbation of the equilibrium distribution f e = f 0e + 5f e , \5f e \ -C fo e , we also 
introduce small perturbations of the electron density and temperature, n e = riQ e + 5n e 
(\Sn e \ <C ^oe), and T e = T 0e + 5T e (\ST e \ T 0e ). In the isothermal approximation, we 
neglect the temperature perturbations, assuming T e = T 0e (with T 0e e^(n e )). Then, 
expanding fo e {n e ,T 0e ) on 5n e and linearizing with respect to the small perturbation, we 
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obtain 



tBGK 



Sfe- 



dfoe(n e ,Ti 



Of: ) 



dn f 



5n P 



with foe(n e ,Toe) denned by ([T]). Differentiating, we obtain 

{2ixhf 



tBGK 



2 5n e e F0 

Vei \ QJe - - — — J0e{n0e,T 0ej 



3 7loe To e 

where e^o = eir(no e ). Finally, noting that 



foe(n 0e , T( 



Oe; 



/oe 
2oe 



1 JOe 



9e ' 



where e = p 2 /2m e is the electron kinetic energy, we obtain 

tBGK 

ei 

with (5n e = J 5f e d 3 p. 
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